Deviation from one-dimensionality in stationary properties and collisional dynamics of 

matter-wave solitons 



Lev Khaykovich 1 and Boris A. Malomed 2 
1 Department of Physics, Bar-Ilan University, Ramat-Gan, 52900 Israel, and 
2 Department of Interdisciplinary Studies, School of Electrical Engineering, 
Faculty of Engineering, Tel Aviv University, Tel Aviv 69978, Israel 

^° . 

By means of analytical and numerical methods, we study how the residual three-dimensionality 
, affects dynamics of solitons in an attractive Bose-Einstein condensate loaded into a cigar-shaped 

C^l ■ trap. Based on an effective ID Gross-Pitaevskii equation that includes an additional quintic self- 

focusing term, generated by the tight transverse confinement, we find a family of exact one-soliton 
solutions and demonstrate stability of the entire family, despite the possibility of collapse in the ID 
equation with the quintic self-focusing nonlinearity. Simulating collisions between two solitons in 
the same setting, we find a critical velocity, V c , below which merger of identical in- phase solitons is 
observed. Dependence of V c on the strength of the transverse confinement and number of atoms in 
the solitons is predicted by means of the perturbation theory and investigated in direct simulations. 
The simulations also demonstrate symmetry breaking in collisions of identical solitons with a nonzero 
phase difference. This effect is qualitatively explained by means of an analytical approximation. 
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I. INTRODUCTION 



It is well known that a trapped atomic Bose-Einstein condensate (BEC) with attractive interactions is stable if the 
number of atoms in it is below a critical value, above which collapse occurs pj. Beneath the collapse threshold, the 
BEC can form stable wave packets in a one-dimensional (ID) "cigar-shaped" trap, which is tightly confined in two 
(transverse) directions, and is unbound along the longitudinal axis. In that case, the stability of bright solitons is 
provided by balance between the quantum pressure, alias matter-wave (MW) dispersion, and mean-field attraction. 
i^i - A single MW soliton and trains of interacting solitons Q have been created in the cigar-shaped optical traps. 
However, while the trapping geometry was nearly one-dimensional, the solitons themselves were far from being ID 
objects. In particular, in Ref. Q{, a stable soliton was only possible if its longitudinal size exceeded the transverse size 
by no more than 20% (note that the situation was affected by an expulsive axial potential, unavoidable in the specific 
experimental setup). Recently, it was shown that the proximity of the soliton to being a 3D object strongly affects 
its properties, such as the character of its motion |4| and interactions lJ-lJ In particular, it was demonstrated that 
a moving soliton immersed in a cloud of thermal atoms is subjected to a temperature-dependent friction force 0. A 
collision between two solitons, which are by themselves stable, in a confined geometry may readily lead to collapse, 
[ if the total number of atoms in the soliton pair exceeds the above-mentioned critical value, and the phase difference 
between the solitons is (close to) zero j(|. The significance of the effective dimensionality of MW solitary pulses is 
further emphasized by the recent observation of formation of a set of nearly-3D mutually repulsive MW solitons (with 
a phase shift of ir between them) as a result of incomplete collapse in an attractive BEC with the number of atoms 
several times larger than the critical value 0,0- 

In addition to being a profoundly important object of fundamental studies, MW solitary waves are also natural 
candidates for applications, such as high-precision atom interferometry and quantum-information processing. Thus, 
a thorough understanding of deviations of their behavior from that of ideal ID solitons is important in this respect 
too. 

In this paper we report results of theoretical investigation of the shape of stationary MW solitons and binary 
collisions between them in the quasi-lD regime, with the aim to identify manifestations of nonsolitonic behavior 
due to the residual multi-dimensionality. The effect of the tightly confined transverse dimensions is taken into 
account through a perturbative self-focusing quintic term added to the corresponding one-dimensional Gross-Pitaevskii 
equation (GPE), as per Refs. and 0. In Section II, we introduce this extended GPE, find a family of its exact 
one-soliton solutions, and demonstrate stability of the entire family, despite the fact that collapse occurs in the ID 
equation with the quintic self-focusing term. In Section III, we investigate soliton collisions within the framework of 
this equation. On the contrary to completely elastic collisions between solitons in the cubic GPE (alias cubic nonlinear 
Schrodinger equation, NLSE), in the presence of the quintic term colliding solitons with zero phase difference, Aip = 0, 
merge into a single pulse if their relative velocity is smaller than a critical value, 2V C . We find the dependence of V c 
on the strength of the transverse confinement and number of atoms in the solitons. For moderate quintic nonlinearity, 
good agreement with an analytic prediction derived from the perturbation theory is found. With a stronger quintic 
term, the numerical results deviate from the perturbation theory, although not dramatically. Finally, we demonstrate 
dynamical symmetry breaking between identical solitons colliding with A(p ^ (in that case, the merger does not 
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occur), as a function of the relative velocity. An explanation to the latter effect is proposed. It is based on estimation 
of a symmetry-breaking parameter, which is a mismatch between the amplitude center and phase center of the soliton 
pair with Aip ^ (exact definitions are given below) . Reasonably good agreement between numerical results and the 
analytical approximation is observed. The paper is concluded by Section IV. 



II. AN EFFECTIVE ONE-DIMENSIONAL GROSS-PITAEVSKII EQUATION AND EXACT SOLITON 

SOLUTIONS 



A. Basic equations 



We start with the standard GPE for a condensate tightly confined in the transverse plane, with the radial coordinate 
r, and unconfined in the axial direction, x: 



ih -£ = -^r v i + Toi UP + o"*>v^ + — — (i) 



where operator acts in the transverse plane, to is the frequency of the trapping potential in this plane, m is the 
atomic mass, and a < is the scattering length. Transition to the quasi-lD description is possible if the change 
of the chemical potential due to the mean-field interaction is much smaller than the level spacing in the transverse 
trapping potential. We briefly recapitulate the corresponding derivation, following, chiefly, Ref. |4|. In the quasi- ID 
limit, the factorized ansatz, ip(r,x,t) — <fi(x,t)x(r,x,t) is used to adiabatically separate fast transverse and slow 
longitudinal dynamics, by neglecting derivatives of x with respect to the slow variables, x and t. By substituting the 
ansatz into Eq. , two decoupled equations are obtained, within the framework of the tight-transverse-confincment 
approximation: 

dq> k 2 d 2 q> 

lh m = -^-¥x- + ^ (2) 



^ n2 1 22 4:irh 2 a l2 , oS 

t l X = — v ±X + —muj r x H n\X\ X, (3) 

2m 2 m 

where the transverse chemical potential, Jx, has to be found from the ground-state solution of Eq. © as a function of 
the ID density, n(x, t) = \4>(x,t)\ 2 . Physical solutions of Eq. (j3Jl exist only if —an < 0.47 0, otherwise transverse 
collapse occurs pj| . In the quasi-lD limit, corresponding to —an << 0.47, the transverse wave function, X: is close 
to the ground state of the 2D harmonic potential, and can be expanded over the set of transverse eigenmodes, </3 m (r): 
x{r,x) = (fo(r) + Tt m C m (x)tp m (r). Coefficients C m are small and can be calculated perturbatively. Accordingly, the 
transverse chemical potential jl can be expanded over powers of the density by means of the perturbative theory, 
jl — fibj + gir>n — §2n 2 + where 

fflD = 2hu;a, g 2 = 24 (in hu;a 2 , (4) 

as shown in Ref. [8j (the subscript "ID" implies that the corresponding coefficient appertains to the standard ID 
model). Substituting the expansion for jl in Eq. one arrives at an effective equation describing the soliton 

dynamics in the quasi- ID limit: 

which is NLSE with the cubic-quintic (CQ) nonlinearity. 

Other approaches to the derivation of the effective ID GPE were also proposed [l2l IT3|. In particular, a more 
complex equation with nonpolynomial (algebraic) nonlinearity was derived, by means of the variational approach to 
the separation of the axial and transverse wave functions, in Ref. |l2j| . Expanding the nonlinearity up to the quintic 
term, one arrives at an equation similar to Eq. J5J, but with a different numerical coefficient. 

NLSEs with the CQ nonlinearity are well known as model equations in nonlinear optics, starting with pioneer works 
Ha. GPEs with the CQ nonlinearity were also used in order to take into account three-body collisions in the BEC 
|17|. However, in the previously considered settings, these equations were always considered with a combination of 
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self-focusing cubic and self-defocusing quintic terms. A drastic difference in the present case is that the quintic term 
is self-focusing [as seen from Eq. @, this conclusion does not depend on the sign of scattering length a, i.e., on 
the self-focusing or defocusing character of the cubic te rm- the same conclusion follows from the expansion of the 
above-mentioned nonpolynomial NLSE derived in Ref. [U- The use of the GPE with the "double-self-focusing" 
CQ nonlinearity, which is the case here, was tacitly assumed impossible, as in this case the equation gives rise to 
collapse. Nevertheless, we will show below that this equation generates meaningful stable solutions. In fact, if the 
cubic nonlinearity is self- focusing, i.e., the scattering length is negative (the case considered throughout the present 
work), the presence of the collapse is a relevant qualitative feature of the effective GPE, as collapse takes place too in 
the full 3D equation, from which Eq. (0 was derived (even if the strong collapse in the full 3D GPE and weak collapse 
in the ID CQ equation bear essential differences). As shown in Refs. |oL ll3L Il5j. the collapse in the 3D equation 
may be avoided under the constraint of N\a\/a± < 0.627, where N is the number of atoms in the condensate, and 
aj_ = \Jnj (mu>) is the harmonic-oscillator length corresponding to the transverse confinement. 



B. Soliton solutions 

Below, we use Eq. @ in the normalized form, 

where gir> < and 32 > are dimensionless interaction constants. In fact, the absolute values of both of them may 
be additionally scaled to be 1, but we find it more convenient to keep these coefficients as free parameters. 

A family of exact soliton solutions to Eq. © can be found as an analytical continuation of the well-known solution 
of the equation with the self-defocusing quintic term ^(| . The result is 



3 ^ 1/4 



( x t) = 2 [ — I e" iM * / - (7) 

\W V vV - Vcosh (2^2jlx) + g ' 



.9 = —x\ —91B , (8) 
2 V 92 

where /i is the soliton's chemical potential that may take any value from < —fi < 00. The squared amplitude of this 
soliton, i.e., the maximum atomic density at its center, is 



A 2 = i,/A( v /7^-.g), 
2 V 92 V ' 



(9) 



and the norm of the soliton, which measures the total number of atoms, is 

f +00 

|2. 




N sol = / \(t>(x)\ 2 dx = W-tan- 1 v ' . (10) 



The soliton's norm and squared amplitude are shown, as functions of the chemical potential, in Fig. ^ 

It is worth to note a drastic difference of this soliton family from its counterpart in the model with the self-defocusing 
quintic term, i.e., gi < 0: in that case, the norm takes all values, < N so \ < 00, while the chemical potential is limited 
to a finite interval, < — fi < |/i| max = 3.9id / 0-® I <72 1 ) - O n the contrary, for the present solution family, Eq. I|10l) 
demonstrates that the norm is limited to a finite interval, 

< N sol < N max = J -?-7r, (11) 

V 852 

while —//is not limited from above (as said above). In fact, iV max in Eq. I|ll|l is a collapse border of the soliton family. 
Further, in the usual CQ model, with g 2 < 0, the amplitude is limited by a finite value, A 2 < 3gm/ (4(? 2 ), while the 

width of the soliton diverges ~ In f (|/x| max — |/i|) -1 ^ at — * |ju| max - I n contrast to this, Eqs. J7J and © show that 

the amplitude of the present soliton family diverges, A 2 s» W— (3/ 32) at // — * —00, and the width of the soliton 
shrinks in the same limit, as l/y/—fj,. This asymptotic behavior of the soliton solution clearly suggests a transition to 
a collapsing solution at jV so i — > A max , see Eq. iJTTjl. 
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FIG. 1: The upper and lower panels display, respectively, the family of exact soliton solutions JjJ, with gm = — 1 and gi = 3/4 
[hence, g — 1, see Eq. (JHJ], in terms of the dependences of the norm and squared amplitude vs the chemical potential, as per 
Eqs. ETUI and 

Equation (|10fl shows that condition dN/dfi < holds for the entire soliton family (see also the upper panel in 
Fig. nj, hence the solitons satisfy the known Vakhitov-Kolokolov (VK) stability criterion [l^. As this criterion is 
only a necessary one, but not sufficient, the stability of the solitons was tested in systematic direct simulations of 
Eq. ©. Results clearly suggest that all the solitons are indeed stable against small perturbations (of course, a large 
perturbation may provoke onset of the collapse). 

It is relevant to mention that the ID GPE with the full algebraic nonlinearity introduced in Ref. ^| gives rise 
to two branches of (implicit) soliton solutions, one stable and one unstable; the branches meet and disappear at the 
point of transition to collapsing solutions. Equation (JHJ) does not give rise to the second branch, as the combination 
of the cubic and quintic terms may be regarded as a truncated expansion of the full algebraic nonlinearity from the 
above-mentioned equation, and this truncation does not pick up the unstable branch. 

It may also be relevant to note that, starting the derivation of the effective ID equation from the 3D GPE with 
the positive scattering length (corresponding to self-repulsive BEC), one will arrive at Eq. JfjJ) with gm > (and 
again with 52 > 0). The corresponding equation, featuring competition between the cubic self- focusing and quintic 
self-defocusing terms, has a family of exact soliton solutions given by the same expressions (|7|)- (|1UII . in which g is 
negative, as per Eq. JSJ. Despite the formal similarity, the latter soliton family is completely different from the one 
presented above. In particular, in the limit of fi — > the solution is not a usual broad small-amplitude soliton, but 
rather an algebraic one, 



The most drastic difference of the soliton family with giD > and <?2 > from the above one is that it features 
dN/dfi > 0, hence this entire family is unstable, according to the VK criterion [algebraic solitons, such as one in Eq. 
I)12|l. are known to be unstable for a different reason 19]]. Besides the fact that all the solitons in the model with 
the positive scattering length are unstable, their physical meaning is doubtful also because the quintic term, which 
appears as a perturbative correction to the cubic one Q, actually dominates over it in these solutions. 




(12) 



III. SOLITON COLLISIONS 



A. Merger of colliding solitons with Aip = 



It is commonly known that collisions between solitons in the one-dimensional NLSE, which is an integrable equation, 
are completely elastic. The force of interaction between the solitons depends on the relative phase between them: 
with Af = and Aip — n, they are attract and repel each other, respectively H3,0]. The quintic term breaks the 
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integrability of the equation, and is expected to make collisions inelastic. For Aip = 0, simulations reveal a critical 
collision velocity, below which two identical solitons merge into a single one. 



1. Analytical considerations 



The merger may be explained by the fact that radiation loss due to the inelastic collision becomes greater than the 
initial kinetic energy of the soliton pair |2l| . This explanation can be implemented in an explicit form if the quintic 
term is treated as a small perturbation. To this end, defining \& = \/\gm\<j>, we rewrite Eq. © in the following form: 

where e = <?2/<7id- I n the zero-order approximation (e = 0), the travcling-soliton solution to Eq. 1|13|) is 

®(x,t) = A sech{A{x - Vt))e- t{ ^- Vx) , (14) 

where A and V are its amplitude and velocity, and the frequency u> = V 2 /2 + \i is a sum of the kinetic energy and 
binding (potential) energy, [i = — A 2 /2, per particle. 

The use of the perturbation theory makes it possible to obtain the following analytical result for the collision 
between solitons with equal amplitudes A, velocities ±V, and a phase shift A(p$ between them [2l|: if the solitons 
are fast, V 2 3> A 2 , the energy loss generated by the emission of radiation during the collision is 

(A£) rad = e 2 A 2 {aA 5 + V 5 e- KV / A [f3 1 cos(A^) + f3 2 sin(A^)]} , (15) 

where a ~ 1381, f3\ ~ 2401 and /?2 ~ 347. Note that the phase-dependent terms are exponentially small. In the same 
approximation, the collision-induced loss of the number of atoms is 

(AAOrad = (21V 2 ) (A£) rad . (16) 

To estimate a merger condition (threshold), we assume that the velocities ±V, which determine the collision-induced 
losses as per Eqs. (|15f) and 11160 . are actually acquired by originally quiescent (or slowly moving) solitons due to their 
mutual attraction (if Atp is close to zero). To this purpose, we note that the effective potential of the interaction 
between far separated identical solitons is, in the case of e = 0, 

U int (X, Aip ) = ~8A 3 e~ AX cos(A(^ ) (17) 

[2fl | , and the effective mass of the soliton is M g = 2 A. In this approximation, the attraction accelerates the two bi- 
phase solitons to self-acquired velocities, ±V se u, that can be found from the energy-balance equation, 2- (M c g V^, lf /2) = 
8A 3 , hence V^cif = 2 A. Substituting this velocity in Eq. l(lT|l shows that the phase-dependent part is less than 10% 
of the phase-independent one, and therefore wc neglect it. Thus, the collision-induced loss of the energy and number 
of atoms (for both solitons) are predicted by the perturbation theory to be 

(A£) rad - ae 2 A 7 , (AiV) rad = (a/2)e 2 A 5 , (18) 

where a is the same numerical coefficient as in Eq. Ijl5|l . 

The energy of a free soliton and its norm (number of atoms), in the e = limit, are 

E so i = ~A 3 + hl cS V 2 , N sol = 2A (19) 

(the negative term in E so \ is the binding energy). First, the norm loss, AN, taken from Eq. (|18|) . gives rise to the 
collision-induced change of the soliton's amplitude: AA = — (AiV) ra d/2 = — (a/4)e 2 A 5 . The corresponding change in 
the binding (potential) energy of both solitons is positive, 

A£ bind = A(~A 3 ) = -2A 2 AA = ^e 2 A 7 . (20) 

Finally, the energy balance predicts a change in the total kinetic energy: 

AE kin = -(A£) rad - AE hind = -(3a/2)e 2 A 7 . (21) 
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FIG. 2: The critical velocity for the merger of colliding solitons, v c /vo (as an experimentally relevant reference value, we take 
vo — 0.21 mm/s), as a function of the strength of the transverse confinement, lo/luo (with ujq — 2tv x 710 Hz, as in Ref. Q)- For 
relatively weak confinement (smaller u>), dependence V c oc lo 2 is observed, as predicted by the perturbation theory (the solid 
line shows the J 2 power law as a guide to the eye) . The dashed line is the power-law fit to the last four points of the numerical 
results, yielding V c oc o; 2 - 29 ± 07 . 



The merger condition states that the loss of the kinetic energy is equal to or exceeds the initial kinetic energy [21j . 
With regard to the expression for the total kinetic energy of both solitons which follows from Eq. I|19|l . = 2AV 2 , 
this condition means that the merger is expected if the initial velocity of each soliton falls below a critical value: 



vr=^A^^K A . .22, 



The derivation of the merger threshold implies that the critical velocity is much smaller than the above-mentioned self- 
acquired velocity, V^ c if = 2v4 (then, the initial velocities of the solitons may be disregarded in the above energy-balance 
analysis, in comparison with V se tf, as it was actually done). Expression indeed satisfies condition V c <C V sc if, as 
e is a small parameter. 



2. Numerical results 



For simulations of soliton collisions in Eq. (jSJ, we chose parameter values close to those in the real experiment 
, where 7 Li atoms were used: a very small scattering length, a — —0.06 nm, transverse oscillation frequency 
lo = 2n x 710 Hz, and the number of atoms N so \ = 4000. However, we did not include any external longitudinal 
potential, in contrast to the expulsive potential that was present in the experiment. Recall that the expulsive potential 
made the soliton stability region very small Q, and actually caused the soliton to be very close to the 3D limit. The 
present simulations do not include the external potential because we are interested not in effects produced by such a 
potential, but rather in small deviations from the one-dimensionality. In fact, a modification of the above-mentioned 
experimental setup, with the aim to make the central segment of the cigar-shaped trap free of any tangible axial 
potential, is quite possible. 

To compare the analytical prediction for the critical velocity, given by Eq. I|22|) . to numerical results, it is necessary to 
express perturbative parameter e in terms of the transverse trapping frequency to. Undoing the above renormalizations, 
one arrives at a conclusion that Eq. (1221) implies a quadratic dependence, V c cx ui 2 , within the framework of the 
perturbation theory. This dependence is indeed observed in simulations at relatively small to, as seen in Fig. J2J. 
However, at larger w, i.e., for stronger transverse confinement, the numerical results feature a greater power in the 
V c (u)) dependence. In particular, the best fit to the last four numerical points in Fig. J3J) yields V c cx C j 2 - 29 ±0- 07 ] which 
demonstrates a small but tangible deviation from the power law corresponding to the perturbation limit. 

We now turn to the dependence of the critical velocity on the number of atoms, N . The analytical prediction, Eq. 
I|22p. clearly implies V c cx iV 3 . In Fig. ©, this dependence is indeed observed at smaller values of N (i.e., for weaker 
nonlinearity) , where the perturbation limit should naturally be valid. Perusal of numerical data shows that, in this 
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FIG. 3: The critical velocity for the merger of two solitons as a function of the number of atoms in each of the colliding solitons, 
iVsoi- At smaller JV so i, i.e., for weaker nonlinearity, the V c oc N^ ol dependence is observed, as predicted by the perturbation 
theory, see Eq. I122B . To say more accurately, the solid curve, which is the fit to the first four numerical points, features a 
power law V c oc N^-? ±0 ' 2 . The dashed curve is the power-law fit to the last four points, showing a different power dependence, 
KociV s 3 ol 62±a05 . 



range, the actual collision-induced radiation loss is very small, and, as a result, the merger does not lead to complete 
fusion of the colliding solitons into a single pulse, but rather to formation of a bound state of two solitons ("weak 
merger"), as can be seen in Fig. 0] Namely, after the first collision, the solitons re-emerge as two distinct wave packets 
which then collide again many times. A similar nearly radiationless inelastic collision, leading to the formation of a 
two-soliton loosely bound state, was recently observed in simulations of a weakly discrete cubic NLSE [22| . 

A definite deviation from the V c oc N^ ol dependence is observed in Fig. [3] for AT so i > 4000, which shows a limitation 
of the perturbative predictions. In this regime of strong nonlinearity, a smooth transition in the collision process 
occurs, from the formation of the above-mentioned long-lived bound state to direct ("strong") merger of two solitons 
into a single pulse, which is accompanied by a burst of radiation. The conspicuous loss of matter with the radiation 
prevents the emerging single pulse from having the number of atoms above the collapse threshold, therefore the pulse 
does not blow up. The transition is expressed in reduction of the life time of the loose bound state before the complete 
merger. In Fig. |SJ which represents the strongest nonlinearity included in the present framework, the bound state 
features only two oscillations. 

For even stronger nonlinearities (which were also considered), the power-law dependence of V c on the number of 
atoms and transverse trapping frequency is observed in the form of iV sol and w 3 , respectively. However, in such an 
extreme regime, the relevance of the quasi-lD model is questionable. In any case, these results convey a clear caveat 
to soliton experiments, in which relative variations in the atom number may be as large as ~ 2: the strong power-law 
dependence of the critical velocity on the number of atoms should be taken into account, to avoid occasional merger 
of solitons. 



B. Symmetry breaking in soliton collisions with Ay> ^ 

We proceed to inelastic collisions of identical solitons with the phase difference of < A(p < ir/2. Numerical 
simulations of Eq. © show a salient effect of symmetry breaking in this case: while the solitons separate after the 
collision, they emerge as two pulses with different amplitudes (then, the velocities are also different, to comply with the 
momentum conservation), as shown in Fig. |SJ It should be mentioned that a similar effect was observed in simulations 
of collisions between identical solitons in some other nonintegrable ID models, chiefly in those describing transmission 
of nonlinear optical pulses, within the framework of the coupled- mode theory, in waveguides equipped with Bragg 
gratings. In that context, the collision-induced symmetry breaking was reported in basic single-core models |23j, and 
in more sophisticated dual-core ones |24j. A similar effect was also observed in collisions between moving solitons in 
the discrete NLSE 0. " 

In order to achieve qualitative understanding of the symmetry breaking, we resort to consideration of an ansatz 
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FIG. 4: Density profiles as a function of time in a regime of "weak" merger (iV so i = 3500, v/vq = 2). After the first collision, 
the two solitons re-appear as two distinct wave pulses which then collide again many times. 




FIG. 5: The same as in Fig. |1] but in the regime of "strong" merger (for iVsoi = 5000, v/vo — 8). A short-lived bound state 
quickly merges into a single breather-like pulse. 



based on a formal linear superposition of two completely overlapping identical solitons J7J) , at some moment of time 
t = to, with velocities ±V and phases ±(1/2)A<^. The ansatz yields the following expression: 



)to 



\ 



^g 2 + A\fi\ cosh (2^%T\x) + g 



cos(Vx + Aip) . 



(23) 



An essential peculiarity of this expression is that the central points of the two last multipliers do not coincide: one 
is found at x — 0, while the other one at x — —Atp/V. This simple observation suggests a concept of the mismatch 
between the amplitude center and phase center of the pair of colliding identical solitons. The mismatch was considered 
as a cause of breaking the symmetry between colliding solitons in the above-mentioned model based on the discrete 
NLSE El. 
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FIG. 6: Density profiles of solitons featuring the symmetry breaking in the collision, for Aip = tv/10. 



To characterize the asymmetry of ansatz <|23[) qualitatively, we introduce its center-of-mass coordinate, 



I-oo a; l ( / , ansatz(^)| 2 rf^ sin(2Ay>) 7 cosh(7u) sinh(7r-u) — tt sinh(7u) cosh(7ru) 
S-™ I ^a„ sa tz (x) 1 2 dx ~ sinh(7rt;) [7 sinh(™) + tt cos(2 Aip) sinh( 7 t;)] ' 



(24) 



with v = V/ and 7 = tan _1 (2-\/|/I|/g). For the qualitative understanding of the situation, we adopt a natural 
conjecture that the strongest possible symmetry breaking is attained at a value of the velocity v = w m ax, which 
corresponds to a maximum of |£| for given Aip. 
and Eq. (|24|l simplifies: 



For the weak quintic nonlinearity <C g 2 ), one has 7 2y^u]/<? 



sin(2A</j) sinh(7ru) — 7i"L> cosh(7ru) 
^/2|/i| sinh(7ru)[sinh(7ru) + 7ru cos(2A<^>)] ' 



(25) 



Asymmetry parameter £ is shown, as a function of V, in Fig. (JJJ by the solid line for A93 = w/10. It characterizes 
the degree of the collision-induced symmetry breaking, and predicts a maximum at some nonzero velocity. Quite a 
similar dependence is indeed produced by numerical simulations of Eq. 10 for the same value of Atp, as shown by dots 
in Fig. Q. The dots display values of the amplitude ratio of the output soliton pair, as found from the simulations. 
Actually, the pulses emerging from the inelastic collisions are breathers with time-dependent amplitudes. Therefore, 
we averaged the amplitudes over long propagation distances after the collision. 

Generally, the numerical data in Fig. JJJ follow the predicted symmetry-breaking parameter. However, two notable 
deviations are observed: a nonsmooth shape of the numerical curve (local modulations near the maximum, which 
tend to give rise to two extra local maxima, definitely exceed an error of the numerical simulations), and a weakly 
decaying tail, which implies that the asymmetry generated by collisions between fast solitons is essentially larger than 
predicted by the above model. Note that the analytical model does not include radiation loss. The latter may enhance 
the asymmetry, as the loss itself is, plausibly, asymmetric too. 

We note that Ref. 0, in which collisions of nearly 3D solitary waves were considered through simulations of the 
full 3D GPE, showed very little symmetry breaking ("population transfer") between colliding solitons with the initial 
phase difference of Aip = tt/10, less than 1% . However, our results predict that the matter transfer (symmetry 
breaking) would be conspicuous at specific values of the collision velocities, which might not be included in the 
analysis reported in Ref. p| 

For very small Aip, we observed chaotic behavior in the output of the collision, similar to what was reported in a 
weakly discrete NLSE [2^ (see also Ref. |2(|). Very recently, chaotic behavior was predicted for collisions of more 
than two MW solitons, in the presence of a longitudinal parabolic trapping potential 27]. In our model, the collision 
between two solitons is sufficient to observe chaotic behavior, which will be reported elsewhere. 
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FIG. 7: Black circles show the numerically found amplitude ratio of two solitons after the symmetry-breaking collision for 
initial phase difference Aip = ir/10. The solid line shows the velocity dependence of the symmetry-breaking parameter, as given 
by Eq. 12511 . for the same value of Aip. 



IV. CONCLUSIONS 

This work aims to understand how the tight confinement in transverse directions affects the longitudinal dynamics 
of matter-wave solitons in the quasi-lD setting. Within the framework of the known model, which reduces the 
multidimensional character of the full Gross-Pitaevskii equation to the appendage of an additional self-focusing quintic 
term to the effective ID equation, we have investigated deviations from the ideal soliton behavior. 

A family of exact stationary solutions for the solitons has been constructed, and it was demonstrated that the entire 
family is stable, despite the possibility of collapse in the modified ID equation (with the negative scattering length). 
We have found inelastic effects in soliton collisions, which are impossible in ideal solitons. Two identical in-phase 
solitons merge into a single pulse, if the collision velocity is smaller than a critical value. In fact, two different types 
of the merger were observed, "strong" and "weak" ones, the former leading to the formation of a loose bound state 
of two solitons that feature repeated collisions, with very weak radiation loss, while the latter means direct fusion 
into a single pulse, which is accompanied by a burst of radiation (in that case, the radiation loss helps the emerging 
pulse to drop the number of atoms below the collapse threshold, and thus avoid the blowup). Both the analytical 
approximation, based on the perturbation theory, and numerical results highlight the strong dependence of the critical 
velocity on the strength of the transverse confinement and the number of atoms in the solitons. Symmetry breaking 
in collision between identical solitons with nonzero phase difference was also found, and partially explain by means of 
the calculation of a phenomenologically denned symmetry-breaking parameter, which measures the mismatch between 
amplitude and phase centers of the colliding solitons. 

This work was supported, in a part, by the Israel Science Foundation, through grant No. 1125/04 (L.K.) and the 
Center-of-Excellence grant No. 8006/03 (B.A.M.). 



[1] J. L. Roberts, N. R. Claussen, S. L. Cornish, E. A. Donley, E. A. Cornell, and C. E. Wieman, Phys. Rev. Lett. 86, 4211 
(2001); J. M. Gerton, D. Strekalov, I. Prodan, and R. G. Hulet, Nature 408, 692 (2000); E. A. Donley, N. R. Claussen, S. 
L. Cornish, J. L. Roberts, E. A. Cornell, and C. E. Wieman, ibid. 412, 295 (2001). 

[2] L. Khaykovich, F. Schreck, G. Ferrari, T. Bourdel, J. Cubizolles, L. D. Carr, Y. Castin, and C. Salomon, Science 256, 
1290 (2002). 

[3] K. E. Strecker, G. B. Partridge, A. G. Truscott and R. G. Hulet, Nature 417, 150 (2002). 
[4] S. Sinha, A. Y. Cherny, D. Kovrizhin, and J. Brand, Phys. Rev. Lett. 96, 030406 (2006). 
[5] N. G. Parker, A. M. Martin, S. L. Cornish and C. S. Adams, |comLma t /0603059 (2006). 

[6] B. B. Baizakov, B. A. Malomed and M. Salerno, Phys. Rev. A 70, 053613 (2004), and an article in Nonlinear Waves: 
Classical and Quantum Aspects, ed. by F. Kh. Abdullaev and V. V. Konotop, pp. 61-80 (Kluwer Academic Publishers: 
Dordrecht, 2004; also available at http://rsphy2.anu.edu.au/~asdl24/Baizakov_2004_61NonlinearWaves.pdf). 



11 



[7] S. L. Cornish, S. T. Thompson and C. E. Wieman, Phys. Rev. Lett. 96, 170401 (2006). 

[8] A. E. Muryshev, G. V. Shlyapnikov, W. Ertmer, K. Sengstock, and M. Lewenstein, Phys. Rev. Lett. 89, 110401 (2002). 
[9] V. M. Perez- Garcia, H. Michinel, and H. Herrero, Phys. Rev. A 57, 3837 (1998). 
[10] M. I. Weinstein, Commun. Math. Phys. 87, 567 (1983); L. Berge, T. J. Alexander, and Y. S. Kivshar, Phys. Rev. A 62, 
023607 (2000). 

[11] L. D. Carr, and J. Brand, Phys. Rev. Lett 92, 040401 (2004); Phys. Rev. A 70, 033607 (2004). 
[12] L. Salasnich, A. Parola, and L. Reatto, Phys. Rev. A 65, 043614 (2002). 
[13] A. Gammal, L. Tomio, and T. Frederico, Phys. Rev. A 66, 043619 (2002). 
[14] Y. B. Band, I. Towers, and B. A. Malomed, Phys. Rev. A 67, 023602 (2003). 
[15] L. D. Carr and Y. Castin, Phys. Rev. A 66, 063602 (2002). 

[16] Kh. I. Pushkarov, D. I. Pushkarov, and I. V. Tomov, Opt. Quant. Electr. 11, 471 (1979); S. Cowan, R. H. Enns, S. S. 

Rangnekar, and S. S. Sanghera Can. J. Phys. 64, 311 (1986). 
[17] Yu. Kagan, A. E. Muryshev, and G. V. Shlyapnikov, Phys. Rev. Lett. 81, 933 (1998); F. Kh. Abdullaev, A. Gammal, L. 

Tomio, and T. Frederico, Phys. Rev. A 63, 043604 (2001). 
[18] M. G. Vakhitov and A. A. Kolokolov, Izv. Vuz. Radiofiz. 16, 1020 (1973) [in Russian; English translation: Sov. J. Radiophys. 

Quantum Electr. 16, 783 (1973)]. 
[19] R. W. Micallef, V. V. Afanasjev, Y. S. Kivshar, and J. D. Love, Phys. Rev. E 54, 2936 (1996). 
[20] V. I. Karpman and V. V. Solov'ev, Physica D 3, 487 (1981). 
[21] Y. S. Kivshar and B. A. Malomed, Rev. Mod. Phys. 61, 763 (1989). 
[22] S. V. Dmitriev and T. Shigenari, Chaos 12, 324 (2002). 

[23] W. C. K. Mak, B. A. Malomed, and P. L. Chu, Phys. Rev. E 68, 026609 (2003); D. R. Neill and J. Atai, Phys. Lett. A 
353, 416 (2006). 

[24] J. Atai and B. A. Malomed, Phys. Rev. E 64, 066617 (2001); Phys. Lett. A 298, 140 (2002); ibid. 342, 404 (2005). 
[25] I. E. Papacharalampous, P. G. Kevrekidis, B. A. Malomed and D. J. Frantzeskakis, Phys. Rev. E 68, 046604 (2003). 
[26] S. V. Dmitriev, P. G. Kevrekidis, B. A. Malomed, and D. J. Frantzeskakis, Phys. Rev. E 68, 056603 (2003). 
[27] A. D. Martin, C. S. Adams, and S. A. Gardiner, cond-mat/0604086 (2006). 



